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ABSTRACT 

We describe the angular power spectrum of unresolved 3.6 ^trnlR light in Spitzer GOODS fields. The amplitude 
of the anisotropy spectrum decreases with decreasing flux threshold to which resolved sources are removed 
from images. When all pixels brighter than a Vega magnitude of 24.6 are removed, the amplitude of the power 
spectrum at arcminute angular scales can be described with an extra component of z > 8 sources with a IRB 
contribution around 0.4 nW m" 2 sr" 1 . The shape of the power spectrum, however, is more consistent with that 
expected for unresolved, faint galaxies at lower redshifts with Vega magnitudes fainter than 23 with a total 3.6 
intensity between 0.1 to 0.8 nW m~ 2 sr" 1 . We confirm this assumption by showing that large-scale power 
decreases rapidly when the unresolved clustering spectrum is measured from a processed HDF-N IRAC image 
where locations of faint ACS sources with no IR counterparts were also masked. Based on resolved counts and 
unresolved fluctuations, we find that, at most, about 7.0 nW m" 2 sr" 1 can be ascribed to galaxies. 

Subject headings: large scale structure of universe — diffuse radiation — infrared: galaxies 



1. INTRODUCTION 

The intensity of the cosmic near-infrared background (IRB) 
is a measure of the total light emitted by stars and galaxies 
in the Universe. While the absolute background has been 
estimated by space-based experiments, such as the Diffuse 
Infrared Background Experiment (DIRBE; Hauser & Dwek 
2001) and the Infra-Red Telescope in Space (IRTS; Mat- 
sumoto et al. 2005), the total measured IRB intensity still 
remains unaccounted for by sources. At 3.6 /im (L-band), the 
total IRB intensity is 12.4 ± 3.2 nW m" 2 sr" 1 (e.g., Wright 
& Reese 2000) and resolved sources in Spitzer lead to about 
5.4 nW m" 2 sr" 1 (Fazio et al. 2004) to ~ 6.0 nW m" 2 sr" 1 
(Sullivan et al. 2006). 

On the other hand, fluctuation analysis of the Spitzer im- 
ages have shown an excess anisotropy, which has been at- 
tributed to first-galaxies containing Population III (Pop-Ill) 
stars at redshifts before reionization (Kashlinsky et al. 2005). 
This possibility is motivated by suggestions in the literature 
that Pop-Ill stars can explain the difference between mea- 
sured and resolved total IRB intensity (e.g., Santos, Bromm 
& Kamionkowski 2002; Salvaterra & Ferrara 2003; Cooray 
& Yoshida 2004). If this were to be the case, then this high- 
redshift component is best studied with intensity anisotropics 
of the IRB since a high-redshift clustering spectrum differs in 
shape from low-redshift clustering (Cooray et al. 2004; Kash- 
linsky et al. 2004). In a recent paper (Sullivan et al. 2006), 
we discussed IR intensity fluctuations from resolved sources 
in Spitzer images down to a Vega mag. of 22.5. In Sullivan 
et al (2006), we also suggested that the anisotropics in the un- 
resolved background (e.g., Kashlinsky et al. 2005) are more 
likely be due to unresolved, faint galaxies at redshifts between 
1 and 4 and not from a new population of high-redshift galax- 
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ies with an IRB intensity of ~ 1 to 2 nW m 2 sr 1 in the L- 
band (Salvaterra et al. 2006; Kashlinsky et al. 2005). 

In this paper, we make a new set of clustering measure- 
ments in the the Great Observatories Origins Deep Survey 
(GOODS; Dickinson et al. 2003). Instead of fluctuations as- 
sociated with resolved light, we measure the anisotropy power 
spectrum in multipole space of the unmasked pixels after re- 
moving resolved sources down to a certain flux limit. With 
pixels removed down a magnitude level fainter than 24.6, we 
find that the background power spectrum is generated by faint 
galaxies at redshifts between 1 to 4, though we cannot estab- 
lish precisely the slope of the faint-end number counts below 
the point source detection level. 

The Letter is organized as following: in the next section, 
we briefly summarize the procedure we used to measure clus- 
tering of unresolved IRB light. In § 3 we present our results 
and discuss implications of our measurements. Throughout 
the paper, we refer to Vega magnitudes (mv eg a = is 280.9 
Jy). For numerical models, we make use of a flat-ACDM cos- 
mology with parameters tt m = 0.3, h = 0.7, and a normaliza- 
tion to the matter power spectrum today at 8 h" 1 Mpc scales 
cr 8 = 0.84. 

2. IMAGING DATA OF GOODS 

Imaging data of GOODS fields (both CDF-S and HDF- 
N) with Spitzer IRAC were first reduced using a point ker- 
nel for drizzling (Fruchter & Hook 1998). In both fields, 
sources were detected and masked out from the images us- 
ing three different techniques. One technique was identifying 
all sources in the SExtractor GOODS IRAC catalog (Bertin 
& Arnouts 1996) that used Mexican hat convolution kernels 
to optimize source deblending and tuned to push the surface 
brightness threshold for detection down to the faintest levels. 
The diameter of the circular mask for each source was scaled 
by the source brightness. Sources which are brighter than 18.3 
mag has a 7.2" diameter mask. This scaled linearly down to 
a 4.8" diameter mask for sources at 19.3 mag. The choice 
of masking diameter is defined by the point spread function 
which has its first Airy bright ring at 2.4" distance from the 
center (the pixel scale is 0.6" in the GOODS IRAC mosaics). 
All sources fainter than 19.2 mag were masked with a 4.8" 
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FIG. 1 . — Original GOODS HDF-N image (i) and various masked versions of the same region (all panels use the same intensity scale). Panels (ii), (iii), (iv), 
and (v) correspond to Images A, B, C, and D, respectively (see Table 1). Panel (vi) is a simulated image of PopIII fluctuations imprinted on image B, which has 
an additional mask to describe locations of faint optical Hubble/ACS sources with no IRAC counterparts. 



Table 1 



Image 


Point Source 


4(7 


Fraction of 


.~SN 
W 




Removal Limit 


flux/njy 


Pixels Used 


(nW 2 m 4 sr 1 ) 


A 


19.2 


105/183 


55% 


(1.52/1.45) x 10"" 


B 


20.2 


72.8/66.4 


46% 


(1.00/1.01) x 10-" 


C 


All 


41.9/41.9 


30% 


(5.97/6.35) x 10~ 12 


D 


B+ACS 


68.6/64.3 


20% 


(7.3/8.1) x 10~ 12 



NOTES. — Parameters of 3.6 /im Spitzer GOODS images used for IRB 
anisotropy measurements. Image C removes all sources down to some 
arbitrary magnitude (see text), while image D includes an additional 
mask for pixels in which faint point sources are present in Hubble/ACS 
GOODS catalogs, but no IR counterparts. The two sets of numbers are 
for HDF-N/CDF-S, and "4 a Limit" is the flux (in njy) of the brightest 
remaining pixel after removing pixels in the image at 4 a above the rms 
measured from unmasked pixel, with clipping only performed once. 



diameter mask. We also generated masks from the SExtrac- 
tor segmentation maps. In this procedure every source has a 
region defined for it in the segmentation map. We grew this 
region by 2 pixels to ensure that the wings of the sources are 
masked as well (see Table 1). In addition, we also generated 
an HDF-N image which masked sources that were present in 
the Hubble/ACS catalogs of GOODS fields (Giavalisco et al. 
2004) but were absent in the IRAC catalog (image D). These 
are typically faint, blue galaxies in B435 and V606 ACS bands 
and must lie at z < 5. Finally, a third procedure was used 
to mask out all sources (image C). This involved smoothing 
the image by a 3x3 boxcar, identifying a sky level by us- 
ing a 15x15 boxcar and finding pixels which are more than 
2 sigma above the sky level, when the width of the aperture 
used to mask out the pixels is 3 pixels wide. 
In each of the images, when measuring clustering, we addi- 



tionally discard pixels with intensities that are 4 a above the 
mean pixel flux level estimated by averaging all background 
pixels that were not masked with the initial source detection. 
This clipping, done once, allows us to create an image where 
the brightest pixel has a flux of 183 njy (image A) to 42 njy 
(image C). Note that given the width of the IRAC PSF, a 
brightest pixel of 42 njy would correspond to a source with 
a Vega magnitude of about 22.1 (0.4 /iJy). The lcr flux un- 
certainty for isolated point sources over most of GOODS 3.6 
/im area is 22 njy and corresponds to a formal signal-to-noise 
of 5 point-source detection limit of 23.5 in Vega magnitudes. 
From simulations, the 50% completeness limit of IRAC Sex- 
tractor catalogs is about 22.5. 

When measuring clustering of the empty pixels or the unre- 
solved IR light, we take Fourier transforms (FT) of the images 
with masked pixels filled with white noise corresponding to 
the rms level of the image. In Thompson et al. (2006), the 
masked pixels were set to an intensity unit of zero. We tested 
both options and found that setting the masked pixels to a zero 
intensity leads to a power spectrum that is larger in amplitude 
by 25% at angular scales corresponding to typical distances 
between large masked areas. This effect is likely negligible 
for the Thompson et al. (2006) analysis, due to the much 
smaller Airy disk of the NICMOS instrument. 

In general, each image is described by an intensity distribu- 
tion at each of the pixels 1(9) where 9 is a two-dimensional 
angular vector on the sky. We take the Fourier Transform of 

the image Ig(£) oc J I(9)e~' e e d 2 9, where I is the wave number. 
The angular power spectrum is constructed as Q = \Ig\ 2 (and 
plotted in Fig. 2). In practice, note that we take discrete trans- 
forms and the procedure is analogous to the process described 
in Appendix A of Thompson et al. (2006). We believe a sim- 
ilar procedure was used in Kashlinsky et al. (2005), but we 
have not been able to implement their exact technique used 
for point source detection and masking. In Fig. 2, we use 
bootstrap sampling to estimate the error bars of the binned 
power spectrum estimates. While we plot the total angular 
power spectrum, at small angular scales the clustering is dom- 
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FIG. 2.— Clustering of the unresolved IRB light in GOODS HDF-N (filled 
symbols) and CDF-S (open symbols for images A and B) fields as a function 
of the magnitude cut in L-band for resolved source detection. For reference, 
we also show the angular power spectrum of anisotropies from Kashlinsky 
et al. (2005) in light cyan color. With decreasing flux cut the amplitude 
of clustering decreases, with the minimum reached in image D. The smaller 
points in green are the cross-correlation between an image made with sources 
in resolved source catalogs with m(3.6fim) > 21 and the HDF-N image C. 
This sets a lower limit on the confusion that arises from residual flux form 
wings of the brighter resolved sources that are not properly removed from the 
images. Finally, smaller points in brown are the cross-correlation between 
C images of HDF-N and CDF-S. This sets the lower limit for any common 
systematics from IRAC detectors. The lines show expected levels of angular 
clustering in the IRB: the top and bottom dot-dashed lines are contributions 
from a population of Pop-Ill containing galaxies at z > 8 with an intensity 
contribution to the IRB of 1.5 nW irT 2 sr -1 and 0.4nW irT 2 sr~' at 3.6 /im, 
respectively (see Cooray et al. 2004 for details). The solid lines show the 
clustering spectrum of resolved and unresolved galaxies based on models for 
number count in Fig. 3 using the conditional luminosity function descriptions 
in Sullivan et al. (2006). The dashed line that scales as I 2 with increasing 
multipole is the measured shot-noise in image C. 



inated by a combination of a finite-density source shot-noise, 
Cf N = J S 2 (dn/dS)dS below the cut-off flux S cut , and instru- 
mental noise. Here, we use the shape of the power-spectrum 
at large I to measure the shot-noise directly from clustering 
spectra and compare with expectations given models of counts 
at the faint-end. 

To address the possibility that the clustering spectra may be 
affected by residual flux from bright sources that are not per- 
fectly masked out, we also cross-correlated images at differ- 
ent magnitude cuts with catalogs of the extracted sources. We 
show the cross-correlation power between detected sources 
fainter than 21st magnitude and image C in Fig. 2. This spec- 
trum is below the power measured in images alone and allows 
us to safely assume that the large-scale clustering pattern is 
not simply a result of residual fluxes next to brighter sources. 
Finally, to estimate any systematics common to the two fields, 
we cross-correlated C images of HDF-N and CDF-S fields 
with imaging arrays centered on a common coordinate sys- 
tem. As shown in Fig. 2, this cross-correlation spectrum is 
well below signal spectra in each of the fields. 

3. RESULTS & DISCUSSION 

In Fig. 2 we compare our measurements with the anisotropy 
spectrum from Kashlinsky et al. (2005). At the faintest mag- 
nitude where sources in the GOODS images are removed (im- 
age C), our fluctuation spectrum agrees with their measure- 




L-band Magnitude (Vega) 

FIG. 3. — (a): Shot-noise level of clustering as a function of source removal 
magnitude (solid lines correspond to counts in panel b) while filled (HDF-N) 
and open (CDF-S) symbols show the measured shot-noise in the 4cr clipped 
images (see Table 1). In addition to the finite density of sources, our measured 
shot-noise also includes a noise floor related to detector noise, which is a 
constant as a function of the magnitude cut. (b): Measured source counts and 
models based on conditional luminosity functions (from Sullivan et al. 2006). 
The shot-noise level in the images with sources masked out agrees with the 
expectation from counts and there is no evidence for a new population of 
sources at magnitudes fainter than point source detection. 



ments, though clustering in GOODS images probe a factor of 
two larger angular scales. At the largest angular scales probed 
by our measurements, the clustering power decreases relative 
to fluctuations presented in Kashlinsky et al. (2005). Such a 
decrease is inconsistent with the description that clustering of 
IR light is associated with z > 8 first galaxies with a 3.6 fj,m 
IRB intensity between 1 to 2 nW m" 2 sr" 1 (dot-dashed line in 
Fig. 2; e.g., Salvaterra et al. 2006). 

To test the origin of fluctuations, beyond the analysis in 
Kashlinsky et al. (2005), we produced images where we addi- 
tionally masked pixels in which faint Hubble/ ACS sources in- 
volving mostly blue galaxies at z < 5 (Giavalisco et al. 2004) 
are present with no IR counterparts. The spectrum of the 
residual background with this image (image D) has the lowest 
amplitude at arcminute angular scales and support the hypoth- 
esis that the excess IR clustering is due to faint galaxies. If we 
simply scale that power spectrum down to match the ampli- 
tude of the clustering spectrum at I around 5000, the implied 
intensity of IR light from such a z > 8 component is about 
0.4 nW m" 2 sr" 1 . There is a large uncertainty in modeling the 
expected clustering of PopIII fluctuations related to gas cool- 
ing and star-formation efficiency among others (see, Cooray 
et al. 2004). To compare with previous interpretations, the 
above intensity is determined with parameters similar to those 
in Salvaterra et al. (2006). 

Since there is a possibility that PopIII fluctuations may be 
removed during the masking and point source removal, we 
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simulated PopIII backgrounds with power spectra given by 
dot-dashed lines in Fig. 2. We added these simulated maps 
(see, Fig. 1 panel vi) to GOODS images and repeated the same 
point source detection and masking procedure The measured 
clustering spectrum of this image recovers the original input 
spectrum of PopIII fluctuations within 10%. The open squares 
in Fig. 2 show the expected spectrum if GOODS HDF-N con- 
tained a PopIII background with an intensity of of 0.4 nW m" 2 
sr" 1 corresponding to the lower dot-dashed line in Fig. 2. At 
multipoles of 10 3 to 10 4 , the power spectrum estimates using 
the combined resolved-source and faint ACS-source mask is 
a factor of 2 higher than measured in the corresponding image 
D. While a careful analysis could further lower the maximum 
intensity of allowed PopIII fluctuations, we can safely state 
that our analysis indicates that the PopIII background in IR 
images at 3.6 pm is well below the previous suggested level 
of 1 to2 nWm" 2 sr" 1 . 

As shown in Fig. 2, the clustering spectra are more con- 
sistent with what is expected for faint galaxies. We model 
expected clustering using the same halo model described in 
Sullivan et al. (2006) for galaxies with ra(3.6/im) > 23. At the 
bright-end, these models are normalized to reproduce 3.6 /im 
luminosity functions from the SWIRE survey (Babbedge et al. 
2006). Since the faint-end shape of the LF is still uncertain, 
that uncertainty leads to a difference in the expected counts 
at the faint end as shown in Fig. 3(b), though even with in- 
creasing counts down to fainter magnitudes, these faint galax- 
ies are still restricted mostly to redshifts below 4 (4 < z < 8 
counts are shown with the green line in Fig. 3b). While the 
measured power spectra are consistent with clustering of these 
faint galaxies, our results do not allow us to conclusively state 
if the number counts flatten or continue to increase with a 
large positive slope at fluxes fainter than point source detec- 
tion. At magnitudes fainter than 23, these counts lead to an 
additional contribution to IR light between 0.1 and 0.8 nW 
m" 2 sr" 1 , for flat and positive-slope counts, respectively. Note 
that the clustering calculations in Fig. 2 are simply an extrap- 
olation from halo models in Sullivan et al. (2006). We have 
not attempted to vary any parameters in our model to obtain a 
better fit. Given the agreement in Fig. 2 between predictions 
and measurements, it is clear that one does not need to invoke 
a new population of high-z galaxies to explain unresolved IR 
light clustering. 

We can also explore the number counts of these galaxies 
with the shot-noise amplitude. In Fig. 3(a) we compare the 
measured shot-noise level from each of the images with the 
same expected under completeness corrected counts shown in 
the right panel. There is a general agreement between ex- 
pected and measured shot-noise level, though again we can- 
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not distinguish between different slopes for counts precisely. 
With improved measurements, it may eventually be possible 
for us to address which is the exact slope of the faint-end num- 
ber counts of IR galaxies. We leave this work for a later study. 
Earlier estimates (Fazio et al. 2004; Sullivan et al. 2006) 
showed that down to a Vega magnitude of about 23, counts 
lead to an intensity of 6.0 nW m~ 2 sr" 1 at 3.6/im. With the 
counts increasing with the positive-slope at the faint end, we 
find that, at most, about 7.0 nW m" 2 sr" 1 of 3.6 /im IRB in- 
tensity comes from known galaxy populations. 

Our results do not completely rule out primordial z > 8 
galaxies with Pop-Ill stars in the IR background. They may 
still be present with an intensity contribution below 0.4 nW 
m" 2 sr" . If they are to explain a larger intensity, then this 
background must be extremely smooth at angular scales be- 
low 10 arcminutes and, for some reason, the clustering spec- 
trum must peak at multipoles lower than 10 3 . Moreover, z > 8 
galaxies are expected to be dominating at lower wavelengths 
than 3.6 /im. Recent clustering measurements by Thompson 
et al. (2006) in the NICMOS Deep Field argue against a large 
contribution to the IR light at 1 .6 /im. In the near future, the 
Cosmic Infrared Background Experiment (CIBER; Bock et 
al. 2006; Cooray et al. 2004) will measure zodiacal light, 
total background intensity, and fluctuations of IR light over 4 
square degrees with rocket-borne imagers and spectrometers. 
It should provide the conclusive answer on whether a residual 
in zodiacal light modeling is responsible for the large unex- 
plained intensity of IR light at these wavelengths. 

To briefly summarize our results and discussion, the clus- 
tering power spectrum of unresolved light in Spitzer GOODS 
images are consistent with faint galaxies below the resolved 
counts and predominantly at redshifts below 5. These faint 
galaxies are expected from extrapolating measured IR LFs, 
but remain unresolved even in deep Spitzer images though, for 
most of them, optical counterparts may exist in Hubble/ACS 
images. They are not a significant contributor to the total IRB 
light, and the 3.6 pm total intensity from resolved sources 
down to a Vega magnitude of 23 and those needed to explain 
unresolved fluctuations is about 7.0 nW m" 2 sr" 1 . 
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